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Abstract. HIPPARCOS trigonometrical parallaxes make 
it possible to compare the location of Tc-rich and Tc-poor 
S stars in the Hertzsprung-Russell (HR) diagram: Tc-rich 
S stars are found to be cooler and intrinsically brighter 
than Tc-poor S stars. 

The comparison with the Geneva evolutionary tracks 
reveals that the line marking the onset of thermal pulses 
on the asymptotic giant branch (AGB) matches well the 
observed limit between Tc-poor and Tc-rich S stars. Tc- 
rich S stars are, as expected, identified with thermally- 
pulsing AGB stars of low and intermediate masses, 
whereas Tc-poor S stars comprise mostly low-mass stars 
(with the exception of 57 Peg) located either on the red 
giant branch or on the early AGB. Like barium stars, Tc- 
poor S stars are known to belong exclusively to binary sys- 
tems, and their location in the HR diagram is consistent 
with the average mass of 1.6 ± 0.2 Mq derived from their 
orbital mass-function distribution (Jorissen et al. 1997, 
A&A, submitted). 

A comparison with the S stars identified in the Mag- 
ellanic Clouds and in the Fornax dwarf elliptical galaxy 
reveals that they have luminosities similar to the galactic 
Tc-rich S stars. However, most of the surveys of S stars 
in the external systems did not reach the lower luminosi- 
ties at which galactic Tc-poor S stars are found. The deep 
Westerlund survey of carbon stars in the SMC uncovered 
a family of faint carbon stars that may be the analogues 
of the low-luminosity, galactic Tc-poor S stars. 
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1. Introduction 

The S stars are late-type giants whose spectra resemble 
those of M giants, with the addition of distinctive molecu- 
lar bands of ZrO (Merrill 1922). Detailed abundance anal- 
yses of S stars (e.g., Smith & Lambert 1990) have shown 
that the overabundance pattern for the elements heavier 
than Fe bears the signature of the s-process nucleosynthe- 
sis (see Kappeler et al. 1989). Furthermore, the C/O ratio 
in S stars is intermediate between that of normal M giants 
(~ 0.2) and that of carbon stars (> 1). S stars were there- 
fore traditionally considered as transition objects between 
normal M giants and carbon stars on the asymptotic giant 
branch (AGB) (e.g., Iben & Renzini 1983). Internal nucle- 
osynthesis processes related to the thermal pulses (a recur- 
rent thermal instability affecting the He-burning shell in 
AGB stars) could possibly produce the s-process elements 
(see the discussion by Sackmann & Boothroyd 1991), and 
the 'third dredge-up' (the descent of the lower boundary 
of the convective envelope into the region formerly pro- 
cessed by the thermal pulse) could bring them to the sur- 
face. The simple M-S-C evolution sequence faces, how- 
ever, several problems, as discussed by e.g. Lloyd Evans 
(1984), Mould & Aaronson (1986) and Willems & de Jong 
(1986). The observation of Tc lines (an element with no 
stable isotopes) in some but not all S stars (Merrill 1952; 
Little et al. 1987) raised another major problem. If the 
s-process really occurred during recent thermal pulses in 
S stars, Mathews et al. (1986) predicted that Tc should 
be observed at the surface along with the other s-process 
elements. 

A breakthrough in our understanding of the evolution- 
ary status of S stars came with the realization that Tc- 
poor S stars are all members of binary systems (Brown 
et al. 1990; Jorissen et al. 1993; Johnson et al. 1993), and 
are likely the cooler analogues of the barium stars, a fam- 
ily of peculiar G-K giants first identified by Bidelman & 
Keenan (1951). Barium stars and S stars share the same 
abundance peculiarities (e.g., McClure 1984) and high rate 
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of binaries with suspected white dwarf (WD) companions 
(Bohm-Vitense et al. 1984; McClure & Woodsworth 1990; 
Jorissen et al. 1997a). Both families are believed to owe 
their chemical peculiarities to the transfer of s-process- 
rich matter from the former AGB star (the progenitor of 
the current WD) to the (main sequence) progenitor of the 
current chemically-peculiar red giant. 

It is therefore currently believed that two very different 
kinds of stars are found among S stars: Tc-rich intrinsic 
S stars that owe their chemical peculiarities to internal 
nucleosynthesis processes occurring during thermal pulses 
on the AGB, and Tc-poor extrinsic S stars, all members 
of binary systems with mass transfer responsible for the 
pollution of the S star envelope. Contrarily to intrinsic S 
stars that need to be thermally-pulsing AGB (TP-AGB) 
stars, extrinsic S stars may populate the first-ascent red 
giant branch (RGB), prior to the core He flash, or the 
early AGB (E-AGB) preceding the TP-AGB, since the 
mass-transfer scenario does not set any constraint on the 
current evolutionary stage of the extrinsic S star. A di- 
rect check of these predictions has been hampered till 
now by the difficulty in evaluating the absolute magnitude 
of S stars. Methods used so far include individual paral- 
laxes (x Cyg: Stein 1991), membership in a binary sys- 
tem with a detected main sequence companion (ir 1 Gru: 
Feast 1953; 57 Peg: Hackos & Peery 1968; T Sgr: Cul- 
ver & Ianna 1975), membership in a moving group (ir 1 
Gru, R Hya: Eggen 1972a; R And, HR 363, o 1 Ori: Eggen 
1972b), membership in a cluster or association (WY Cas: 
Mavridis 1960; TT9, TT12: Feast et al. 1976), Call K-line 
emission width (HD 191630, 57 Peg: Warner 1965; o 1 Ori: 
Boesgaard 1969). Most of these individual estimates were 
used by Scalo (1976) to locate S stars in the Hertzsprung- 
Russell (HR) diagram, with the conclusion that they fall 
above the luminosity threshold for TP-AGB stars, as de- 
rived from the early AGB model calculations available to 
Scalo. Although Scalo's conclusion appears to support the 
customary M-S-C sequence, statistical estimates based on 
the kinematics or space distribution of S stars pointed out 
a systematic difference between the average luminosity of 
non- or weakly-variable S stars and of Mira S stars, the 
latter being about 3 visual magnitudes brighter than the 
former (Takayanagi 1960; Yorka & Wing 1979). Such a 
segregation apparent in these early studies might already 
hint at the current dichotomy between extrinsic and in- 
trinsic S stars. 

The rather heterogeneous set of former luminosity 
determinations for S stars, along with our revised un- 
derstanding of their evolutionary status, prompted the 
present study, based on the trigonometrical parallaxes pro- 
vided by the HIPPARCOS satellite. Its aim is to locate 
both kinds of S stars in the HR diagram, and to compare 
their respective locations with stellar evolutionary tracks. 

Our HIPPARCOS sample of S stars is described in 
Sect. ||. The method used for deriving their bolometric 
magnitudes is detailed in Sect. ||. The HR diagram of S 



stars is discussed in Sect. || with special emphasis on the 
comparison with theoretical evolutionary tracks and with 
S and C stars in the Magellanic Clouds. Finally, the cor- 
relation between the infrared colours of S stars and their 
location on the giant branches is discussed in Sect. ||. 

2. S stars in the HIPPARCOS catalogue 

A full discussion of the HIPPARCOS mission can be 
found in the HIPPARCOS catalogue (ESA 1997). Cross- 
identifications between the HIPPARCOS catalogue and 
the General Catalogue of Galactic S Stars (Stephenson 
1984; GCGSS) yield 63 stars in common between the two 
catalogues. Of these, 22 are part of the HIPPARCOS gen- 
eral survey (a systematic monitoring of all stars down to 
magnitudes 7.3 - 9, depending on galactic latitude and 
spectral type), and 41 belong to samples included in the 
HIPPARCOS Input Catalogue (Turon et al. 1992abc) for 
a particular purpose. The sample of S stars studied in 
this paper can thus in no way be considered as a complete 
sample, since it rather reflects the particular interests pre- 
vailing at the time of construction of the HIPPARCOS In- 
put Catalogue. A comparison with the GCGSS indicates 
that the HIPPARCOS sample of S stars is nevertheless 
complete down to V — 7.5. 

Table lists various identifications of the S stars con- 
sidered in this paper. The different columns contain the 
following data: 

1. HIPPARCOS Input Catalogue (HIC) number; 
2-3. GCGSS and HD numbers; 

4. variable name in the the General Catalogue of Variable 
Stars (Kholopov et al. 1985; GCVS); 

5. variability type from the GCVS; 

6. spectral type from the GCGSS; 

7. presence (y) or absence (n) of technetium lines (from 
Little et al. 1987, Smith & Lambert 1988, from the compi- 
lation of Jorissen et al. 1993, or from Van Eck & Jorissen 
(in preparation) in the case of HIC 28297); 

8-9. annual parallax ir and its standard error e>v (both 
expressed in milliarcseconds) . 

Table || lists the basic photometric properties of the 
S stars in the HIPPARCOS sample. The various columns 
contain the following data: 
1. HIC number; 

2-3. median Hp magnitude and its standard error UHp; 

4. colour excess £"b-v; 

5. Johnson Vj magnitude from the HIPPARCOS cata- 
logue; 

6. dereddened (V — K) colour index; 

7. dereddened Kq magnitude; 

8. bolometric correction BCk to the K magnitude; 

9. bolometric magnitude Mboi (see Sect. || for details 
about the derivation of the quantities listed in columns 
4-9); 

10. method used for deriving the bolometric correction: 
(1) by integrating under the energy curve, or (2) from the 
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{BC K , V-K) relation of Bessell & Wood (1984); 

11. source for the HIPPARCOS (V - J)c colour [A-G 
from VRIc photometry; H-K: from UBVj photometry 
L-P: from HIPPARCOS and Star Mapper photometry 
Q: specific treatment applied to long-period variables; R 
from spectral type; see the HIPPARCOS and Tycho Cat- 
alogues (ESA 1997, Vol. 1, Introduction and Guide to the 
Data) for details]; 

12. duplicity flag (field H59 of the HIPPARCOS catalogue; 



see Sect. |4T 

13. other identifications ('Hen' refers to the survey of S 
stars by Henize 1960). 
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Fig. 1. The spectrum of HIC 42650 = GCGSS 544, with the 
principal spectral features identified 



One star (HIC 42650 = GCGSS 544) appears outlying 
in many respects. With the largest parallax in the sample 
and the second faintest Hp magnitude, it is intrinsically 
much fainter than all the other stars in the sample. It is 
also somewhat bluer, with (B — V)j = 1.39 and (V—I)c = 
1.65, compared to average colour indices of {(B — V)j) = 
1.79 and {(V - I)c) = 2.72 for the whole sample. A low- 
resolution spectrum (AA ~ 0.3 nm) of this star, covering 
the spectral range 440 - 820 nm, has been obtained at the 
European Southern Observatory (ESO, La Silla, Chile) on 
the 1.52-m telescope equipped with the Boiler & Chivens 
spectrograph (grating #23 + filter GG 420; 114 A mm -1 ) 
and a Loral/Lesser thinned, UV-flooded 2048x2048 CCD 
(CCD#39; 15 /*m pixels). The flux response curve of the 
system has been calibrated using the spectrophotometric 
standard star LTT 4816. 



The spectrum of GCGSS 544 appears to be that of a 
M0-M1 dwarf, as can be seen on Fig. [TJ from the strong 
MgH (A478.0 and A521.1) and CaH (A638.2 and A638.9) 
bands, as well as from the strong Na D line superimposed 
on moderately strong TiO bands (Jaschek & Jaschek 
1987). This classification is consistent with its absolute 
visual magnitude M v = 8.6 derived from its HIPPARCOS 
parallax and Vj magnitude. That star is thus misclassificd 
as S star, since there is no trace whatsoever of ZrO bands 
in its spectrum. 

The case of T Cct (= HIC 1728) is also conflicting, 
since it was classified as M5-6Se in the original paper by 
Keenan (1954) defining the S class, and reclassified as M5- 
6Ib-II in the Michigan Spectral Survey (Houk & Cowley 
1975). The spectra of supergiants and weak S-type stars 
look very similar at the low plate dispersions used in classi- 
fication work, and are therefore easily confused (e.g., Lloyd 
Evans & Catchpole 1989). T Cet has nevertheless been 
kept in our final list until higher resolution spectra resolve 
these equivocal classifications. 

The star HIC 27066 (= GCGSS 157) has HIPPARCOS 
colours that do not match those of an S star \{B — V)j = 
0.80 and (V - I)c = 0.83]. HIPPARCOS found it to be a 
close visual binary with a separation of 0.262 arcsec, the 
companion being 1.55 mag fainter in the Hp band. The 
colour indices measured by HIPPARCOS correspond to 
the composite light from the system. 

The only other S star found by HIPPARCOS to be 
a close visual binary (with an angular separation of 0.18 
arcsec) is HIC 19853. Although its colour indices are not 
atypical for a late-type star, the measured colours must 
clearly be composite since the companion is only 0.35 mag 
fainter in the Hp band. 

These two stars have been excluded from our final sam- 
ple. Among the 61 remaining S stars, 20 do exhibit tech- 
netium lines in their spectrum, 21 do not, and the Tc 
content is unknown for the remaining 20 stars. 

The only trigonometrical parallax for an S star with a 
small relative error available in the literature prior to the 
HIPPARCOS mission is that obtained by Stein (1991) for 
the Mira S star \ Cyg- His value (8.8 ±1.9 mas) is in 
good agreement with the HIPPARCOS parallax listed in 
Table % 

3. Colours and bolometric magnitudes 

3.1. Colours 

Although the (B — V)j and (V— I)c (where the subscripts 
J and C stand for the Johnson and Cousins photometric 
systems, respectively) colour indices are directly provided 
by the HIPPARCOS catalogue, we felt that the (V - K)j 
colour index, derived on a case-by-case basis, is a more 
appropriate temperature indicator. For most of the stars 
in our sample, the (V — I)c index listed in the HIPPAR- 
COS catalogue has been derived from a fiducial relation 
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Table 1. S stars in the HIPPARCOS catalogue: Identifications and parallaxes 



HIC 


GCGSS 


HD 


Var 


Var type 


Sp. 


Tc 


TV 




621 


3 


310 






S3,l 


n 


2.50 


0.69 


1728 


8 


1760 


T Cet 


SRc 


M5-6Se; M5-6Ib-II 


y 


4.21 


0.84 


1901 


9 


1967 


R And 


M 


S5-7/4-5e 


y 


-0.06 


6.49 


5091 


22 


6409 






M2wkS 


n 


2.43 


0.89 


5772 


26 


7351 






S3/2 


n 


3.21 


0.82 


8876 


45 








S3/1 


n 


-1.92 


1.50 


10687 


49 


14028 


W And 


M 


S7/le 


y 


-1.17 


3.17 


17296 


79 


22649 


BD Cam 


Lb 


S4/2 


n 


6.27 


0.63 


21688 


104 


29704 






S: 


n 


1.84 


1.00 


22667 


114 


30959 


o 1 Ori 


SRb 


S3/1 


v 

j 


6.02 


0.94 


25092 


133 


35155 






S4,l 


n 


1.32 


0.99 


26718 


149 


37536 


NO Aur 


Lc 


M2S 


y 


2.38 


0.97 


28297 


178 


40706 






S2,l 


n 


1.17 


0.88 


30301 


212 


44544 


FU Mon 


SR 


S7/7 (SC) 




0.30 


1.58 


32627 


260 


49368 


V613 Mon 


SRb: 


S3/2 


n 


1.65 


1.11 


32671 


265 


49683 






M4S 




-0.18 


1.07 


33824 


283 


51610 


R Lyn 


M 


S5/5e 
S4*l 




-3.39 


1.79 


34356 


307 


53791 


R Gem 


M 


v 
y 


-6.22 


6.50 


35045 


312 


54587 


A A Cam 


Lb 


M5S 


V 


1.24 


1.02 


36288 


347 


58521 


Y Lyn 


SRc 


M6S 


v 

J 


4.03 


1.33 


37521 


382 


61913 


NZ Gem 


SR 


M3S 


n? 


3.19 


0.79 


38217 


411 


63733 






S4/3 


y ■ 


0.00 


0.99 


38502 


422 


64332 


NQ Pup 


Lb 


S5/2 


y 


3.01 


1.11 


38772 


436 




SU Pup 


M 


S4,2 




-1.75 


1.81 


40Q77 


AQA 


7D97R 


V V_/I1C 


M 


DO/ UU 




9fi ^8 


49 7 A 

4Z. ( 4 


4:<JU<JO 




7871 9 


1XO V_/I1C 




MRS 


y 


8.21 


Q8 




799 
i zz 








M3 


n 


2.10 


u.yu 


CQQ/I A 

oyo44 


788 




1511 oru 








1 PA 
1.04 


n qq 


OZ1ZO 


ono 


i i no i o 
110813 


b UMa 


1V1 


bo/be 


y 


rt c?o 
0.b3 


n n a 
0.94 


a a a i o 


olo 








coo 
BO, 3 




i nn 
-1.90 


Z.99 


64778 


816 


115236 


U Y Oen 


CD 

SK, 


So/8 (OS) 




1.66 


1.04 


00 ( oo 


80£ 
OZO 


lloDoO 






oo,z 


n 


0.4Z 


U.OD 


1 U 1 U 


890. 

ozy 








IVllWKD 


n 


9 97 
Z.Z / 


1 (\p 


68837 






U Cir 


SR 


c 




-1.92 


1.51 


71348 


804 


110994 


BQ Oct 


Lb: 


S5 1 




2.08 


0.57 


79Q&Q 


8IS7 


131217 




ms'9 


n 


4.24 


1.02 




yuo 


142143 


ST Her 


SRh 


MP, ^ 

1V1U.OO 


y 


3.22 


n 7^ 


SI Q70 










M2S 




^ ^9 
o.oz 


1.14 




938 


151011 






Swk 




4.30 


1.07 


87850 




1 63990 


OP Her 


SRb 


M6S 


y 


3.26 


0.54 


88620 


1014 


164392 










2.09 


1.06 


88940 


1023 


165774 






S4,6 


n 


0.23 


1.50 


89316 


1025 


165843 






S2,l 




1.47 


1.05 


90723 


1053 


170970 






S3/1 


y 


1.83 


0.67 


94706 


1117 


180196 


T Sgr 


M 


S5/6e 


y 


-31.67 


9.28 


97629 


1165 


187796 


X Cyg 


M 


S7/1.5e 
S6/3 


y 


9.43 


1.36 


98856 


1188 


190629 


AA Cyg 


SRb 


v 


0.86 


0.88 


99124 


1192 


191226 




MlS-M3SIIIa 


n 


0.39 


0.71 


99312 


1194 


191589 






S: 


n 


2.25 


0.77 


99758 


1195 


191630 






S4,4 


y 


1.18 


0.81 


100599 


1211 




V865 Aql 


M 


S7,2 




-1.28 


1.91 


101270 


1224 


195665 


AD Cyg 


Lb 


S5/5 




-0.96 


1.15 


103476 


1254 


199799 




MS 




2.16 


0.82 


110146 


1292 


211610 


X Aqr 


M 


S6,3e: 




-4.01 


5.73 


110478 


1294 


212087 


tt 1 Gru 


SRb 


S5,7: 


y 


6.54 


1.01 


112227 


1304 


215336 






Swk 


11 


0.74 


0.92 


112784 


1309 




SX Peg 


M 


S3/6e 




2.12 


2.99 


113131 


1315 


216672 


HR Peg 


SRb 


S4/1 


y 


3.37 


0.94 


114347 


1322 


218634 


57 Peg 


SRa 


M4S 


n 


4.28 


0.88 


115965 


1334 






S2/3: 


n 


1.72 


1.26 


HIPPARCOS close visual binaries 












19853 


89 


26816 






S 


y 


3.79 


1.05 


27066 


157 








S 




1.13 


3.50 


Misclassified S star 
















42650 


544 








MS a 




32.13 


1.50 



Comparison barium star 
68023 121447 K7IIIBa5; SO n b 428±71 pc c 



Remarks: 

a: HIC 42650 is an early M dwarf rather than a MS star 
b: Little et al. (1987) 
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Table 2. S stars in the HIPPARCOS catalogue: Photometry 



HIC 


Hp 


0"Hp 


E B -v 




V - K) 


K 


BC K 




i 


Sv-i 


Dupl 


Other ident. 


621 


7.570 


0.003 


0.000 


7.50 


4.38 


3.12 


2.74 


-2.14 


2 


O 




Hen 1 


1728 


5.439 


0.033 


0.000 


5.61 


6.5 


-0.89 


2.91 


-4.85 


1 


o 






1901 


10.705 


0.011 




10.71 




-0.11 


3.49 




2 


K 


V 




5091 


7.462 


0.004 


0.042 


7.44 


5.04 


2.27 


2.84 


-2.95 


2 









5772 


6.395 


0.004 


0.049 


6.33 


4.65 


1.53 


2.68 


-3.25 


1 


F 




HR 363 


8876 


9.072 


0.002 




8.97 










2 


L 




+2T255 


10687 


8.075 


0.226 




8.60 




0.7 


3.49 




2 


O 






17296 


5.095 


0.003 


0.041 


5.06 


4.78 


0.15 


2.77 


-3.09 


1 


C 


o 


HR 1105 


21688 


8.225 


0.009 


0.000 


8.23 


5.64 


2.59 


2.91 


-3.16 


2 


o 




Hen 3 


22667 


4.702 


0.005 


0.121 


4.71 


4.96 


-0.62 


2.85 


-3.87 


1 


o 


V 




25092 


6.882 


0.011 


0.128 


6.82 


4.46 


1.97 


2.76 


-4.66 


2 


F 






26718 


6.243 


0.011 


0.245 


6.23 


4.67 


0.8 


2.74 


-4.57 


1 


c 






28297 


8.987 


0.002 


0.029 


8.88 










2 


L 




Hen 7 


30301 


8.520 


0.032 


0.327 


9.76 


7.35 


1.4 


3.03 




1 









32627 


7.749 


0.004 


0.051 


7.73 


4.93 


2.64 


2.83 


-3.43 


2 


O 






32671 


8.232 


0.006 




8.38 




2.07 


3.49 




2 


O 






33824 


9.922 


0.029 




9.93 




1.91 


3.49 




2 


K 






34356 


7.529 


0.456 




7.53 




2.12 


3.49 




2 


P 


V 




35045 


7.578 


0.007 


0.039 


7.69 


6.01 


1.56 


2.95 


-5.01 


2 


O 


V 




o r* c\ c\ o 

36288 


6.897 


0.005 


0.058 


7.27 


7.63 


-0.54 


3.08 


-4.43 


1 


O 


V 




37521 


5.587 


0.003 


0.027 


5.55 


4.85 


0.61 


2.86 


-4.01 


1 






TTT) ori/ > '7 

HR 29o7 


OOOI 7 

oozl ( 




a aao 
U.UUz 




7 AA 

/ .90 






o ao 
3.4y 




o 
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oooUz 
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(V" — I)c= f(B — V)j calibrated on normal M giants (as 
indicated by the flags H-P in column Sy-i of Table |J). 
Several S stars in our sample, however, fall outside the 
validity range of such a calibration, so that (V — I)c as 
provided by the HIPPARCOS catalogue is not a reliable 
temperature indicator for our purpose. 

As discussed by Ridgway et al. (1980), B — V for late- 
type stars is not a good temperature indicator either, 
because the temperature has opposite effects on B — V: 
when the temperature decreases, the absorption in the V 
band due to TiO bands increases (thus decreasing B — V), 
whereas the black body continuum tends to increase the 
B — V index. By contrast, these two effects act together 
to increase V — K. 

The V — K index has therefore been constructed 
from individual K magnitudes collected from the liter- 
ature, mostly from the Two-Micron Sky Survey (Neuge- 
bauer & Leighton 1969; TMSS), but also from Wing & 
Yorka (1977), Catchpole et al. (1979), Chen et al. (1988), 
Noguchi et al. (1991) and Fouque et al. (1992); it is avail- 
able for 46 stars as listed in Table ||. The adopted V mag- 
nitude corresponds to the Vj magnitude listed in the HIP- 
PARCOS catalogue. It has been obtained by the HIPPAR- 
COS reduction consortia from a (H p — Vj, (V — I)c) rela- 
tion, with (V — I)c derived by different methods as listed 
in column S Y -i of Table [§ According to the HIPPARCOS 
reduction consortia (ESA 1997, Vol. 1, Introduction and 
Guide to the Data) , that relation is quite well defined from 
classical photometry in the range —0.4 < (V — I)c< 3.0 
(which holds for 47 S stars), with uncertainties of less than 
0.01 mag. The red extension down to (V — I)c= 5.4 (15 
S stars) was defined using observations of Mira variables 
devoted to that purpose, resulting in an uncertainty of 
0.03-0.05 mag. For indices down to (V — I)c= 9.0 [only \ 
Cyg is concerned, with (V - 7) c = 6.1], the HIPPARCOS 
reduction consortia adopted a linear extrapolation of the 
previous relation. To avoid the uncertainties inherent to 
such an extrapolation, a (flux) average of the numerous 
Vj measurements of x Cyg available in the literature has 
been preferred over the linear extrapolation. 

In summary, the uncertainty on Vj introduced by the 
colour transformation applied on Hp is largely offset by 
the fact that the Hp magnitude is a good time average 
over a uniform period of time (the duration of the HIP- 
PARCOS mission), the same for all the stars. As discussed 
in Sect. 



4.1 



the intrinsic variability of S stars is indeed a 
major source of uncertainty on their location in the HR 
diagram. 

3.2. Bolometric magnitudes 

Since most of the flux from S stars is radiated in the near- 
infrared, the bolometric correction is best determined from 
near-infrared magnitudes. Moreover, the near-infrared and 
bolometric variability is much smaller than the visual am- 
plitude of variations (Mira-type variables which vary typ- 



ically by 6 to 8 mag in the visual region, vary by less 
than 1 mag in the K band; Feast et al. 1982). Therefore, 
the bolometric correction BCk in the K band (defined as 
M\y Q \ = K + BCk) has been adopted in this work (and 
listed in Table |). 

In order to compute bolometric corrections, an exten- 
sive set of magnitudes ranging from the ultraviolet to the 
far IR has been collected from the literature. This set in- 
cludes Johnson U BV RIJ H K LM magnitudes when avail- 
able, as well as good quality fluxes at 4.2, 11.0, 19.8 and 
27.4 ^m from the Revised Air Force Four-Color Infrared 
Sky Survey (Price & Murdock 1983), and four-colour in- 
frared photometry from Gillett & Merrill (1971). The 
IRAS 12, 25, 60 and 100 ^m fluxes (with a quality flag 
3) from the second edition of the Point Source Catalogue 
(IRAS Science Team 1988) were also used, or when avail- 
able, the reprocessed IRAS fluxes provided by Jorissen & 
Knapp (1997). When several measurements in the same 
filter are available, a flux average has been computed. 

All magnitudes bluer than 4.8/xm have been corrected 
for interstellar reddening and absorption, using the ex- 
tinction law as provided by Cohen et al. (1981) for the 
BVRIHKL filters and by Koornneef (1983) for the J 
and M filters. The visual extinction Ay was derived ei- 
ther from Neckel & Klare (1980) with the distance derived 
from the HIPPARCOS parallax, or from Burstein & Heiles 
(1982) for stars with galactic latitudes \b\ > 10°. In the lat- 
ter case, the -Eb-v value provided by Burstein & Heiles 
(1982) was reduced by the factor [1 — exp(— 10c?sin |6|)], 
where d stands for the distance in kpc. In the re- 
maining cases, the cosecant formula (Feast et al. 
1990) £b-v = 0.032 (cosec|6| - 1) [1 - exp(-10 dsin |6Q1 
was used. The adopted Eb-v values are listed in Table |[ 

The deredenned BV RI J K LM magnitudes have then 
been converted into fluxes using the zero-magnitude fluxes 
provided by Johnson (1966). This particular choice will be 
justified below. The zero-magnitude flux in the H band 
was taken from Jaschek (1978) and, for the remaining 
bands, from the original papers quoted above. 

A limited sample of 17 stars have enough broad-band 
colours available (BVRIJHK, as well as IRAS 12, 25 
and 60 ^m) to derive the bolometric magnitude by a di- 
rect integration of the available fluxes over wavelength. 
More precisely, the trapezoidal rule has been used on the 
curve XF\ versus log A. The zero point of the bolomet- 
ric magnitude has been defined from the requirement that 
L = 3.86 10 33 erg s _1 corresponds to Afboi = 4.75 for the 
Sun. 

If the available photometric data was too scarce to 
derive the bolometric magnitude from a direct integration, 
it has been derived instead from the (BCk, V—K) relation 
of Bessell & Wood (1984) applicable to oxygen-rich stars. 
For the 17 S stars where both methods are applicable, 
they yield consistent results (with a r.m.s. deviation of 0.1 
mag) , provided that the zero-magnitude fluxes of Johnson 
(1966) be adopted (as was done by Bessell & Wood 1984). 



S. Van Eck et al.: The HR diagram of S stars 



7 



If the zero-magnitude fluxes listed by Jaschek (1978) are 
used instead, somewhat lower bolometric corrections BCk 
(i.e. brighter bolometric magnitudes, by about 0.06 mag) 
are obtained. 

A bolometric magnitude Mb Q ] = —3.74, based on the 
same HIPPARCOS parallax, has recently been derived by 
van Leeuwen et al. (1997) for x Cyg, and agrees well with 
our value —3.58 listed in Table ||. 

4. The HR diagram of S stars 

The V — K colours and the bolometric magnitudes derived 
as discussed in Sect. 0, combined with the HIPPARCOS 
parallaxes, provide the HR diagram of S stars presented in 
Fig. ^. Only the 30 S stars with an available K magnitude 
and with < a^/n < 0.85 (see Fig. ||) have been plotted. 
Despite the sometimes large un cert ainties affecting V — K 
or Mboi (as discussed in Sect, fhl] below), a segregation 
between extrinsic and intrinsic S stars is readily apparent, 
with extrinsic S stars being intrinsically fainter and bluer 
than intrinsic S stars. 



-5 



-4 
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Ah H 



^ Tc— rich 

□ Tc-poor 

• Tc unknown 

* HD 121447 



-1 



8 



(V-K) 



Fig. 2. The HR diagram for S stars with < a n /ir < 0.85. 
Filled triangles correspond to Tc-rich S stars, open squares 
to Tc-poor S stars and dots to S stars with unknown Tc. 
HD 1214 47, the boundary case between Ba and S stars (see 
Sect. 4.2), is represented by *. The error bar provides the un- 
certainty on V — K caused by the intrinsic variability of the 
HIPPARCOS Hp magnitude; it covers the range in Hp be- 
tween the 5th and 95th percentiles 




Fig. 3. Distribution of the relative error a-nfir on HIPPARCOS 
parallaxes n for S stars. The hatched and black parts of the his- 
togram correspond to Tc-poor and Tc-rich S stars, respectively. 
Stars with negative parallaxes and stars with a-n/n > 0.85 have 
been assigned to the leftmost and rightmost bins, respectively 



4-1. Uncertainties on the HR diagram 

The two major sources of uncertainty on the position of 
an S star in the HR diagram are the stellar intrinsic vari- 
ability and the error on the parallax it. 

The intrinsic variability of S stars has an impact on 
both V — K and Mboi- The variability in the V band and, 
to a lesser extent, in the K band, leads to a variation of 
V - K, and thus of BC K (Sect. If BC K was derived 
from simultaneous V and K measurements, it could be ex- 
pected that the variations in K and in BCk would cancel 
to a large extent, thus leaving only a moderate variation 
in Mboi- However, since the V and K data used in this 
study do not result from simultaneous observations, the 
expected cancellation will not occur. Its impact on Mboi 
is, however, difficult to evaluate. It seems nevertheless un- 
likely that this effect could lead to a systematic upwards 
shift of all intrinsic S stars that would be responsible for 
the observed segregation between intrinsic and extrinsic S 
stars. The impact of the intrinsic variability on V — K is 
easier to estimate. Its effect is shown on Fig. || from the 
variation recorded in the Hp magnitude over the ~ 1230 d 
duration of the HIPPARCOS mission. The smaller varia- 
tion due to K has not been taken into account. Clearly, 
the uncertainty on V — K due to the intrinsic variabil- 
ity does not jeopardize the observed segregation between 
extrinsic and intrinsic S stars. An interesting by-product 
of Fig. |^ is the increase of the amplitude of variations to- 
wards cooler and more luminous stars (see also Eyer &: 
Grenon 1997 and Jorissen et al. 1997b). 

The impact on the bolometric magnitude of the un- 
certainty on the parallax is shown in Fig. ^, presenting 
(separately for intrinsic and extrinsic S stars) the range of 
Mboi corresponding to n±a 7T (see Arenou et al. 1995 and 
Lindegren 1995 for a discussion of the external errors of 
HIPPARCOS parallaxes and of the accuracy of the zero- 
point). As discussed by van Leeuwen et al. (1997), there is 
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Fig. 4. The uncertainty on Mboi due to the uncertainty ov 

on the parallax n. The error bar extends from Mboi(7r + ov) to 

Mboi(7r — ov). Left panel: parallactic errors for Tc-poor S stars; 

right panel: parallactic errors for Tc-rich S stars. Symbols are 

as in Fig. ^ 



a specific error source on the parallax for nearby Mira vari- 
ables, as some of these stars were found to have asymmet- 
rical spatial light distributions. Changes in these asym- 
metries might affect the derived parallax. The flag 'V in 
field H59 of the H1PPARCOS catalogue possibly reflects 
such effects, as it refers to a 'variability-induced mover'. 
This flag is set for 14 stars in our sample, mostly nearby 
Miras (see column 'DupP in Table |^). Like van Leeuwen 
et al. (1997), we assume that any such effect will add only 
additional random scatter to the mean results. 

Part of the overlap in luminosity between extrinsic and 
intrinsic S stars in the HR diagram may actually be at- 
tributed to the large error bars of the interloping stars 
(NQ Pup and HD 170970o among intrinsic S stars, and 
HD 35155 among the extrinsic S stars). That explana- 
tion does not hold true, however, for the high-luminosity, 
Tc-poor S star 57 Peg (Mboi = —4.3), which has a small 
uncertainty on its parallax (p^jix = 0.21). That star is 
special in many respects, since it has an A6V compan- 
ion instead of the WD companion expected for extrinsic S 
stars in the framework of the binary paradigm (Sect. 1). 
It is further discussed in Appendix A. 

The case of HD 35155 deserves further comments. This 
extrinsic S star has the second largest relative error on 
the parallax. It is a binary system with an orbital period 
of 642 ± 3 d in a nearly circular orbit (Jorissen & Mayor 
1992), yielding an orbital separation of about 2 AU assum- 
ing typical masses of 1.6 and 0.6 M Q for the S star and 
its suspected white dwarf companion (see Jorissen et al. 
1997a). The corresponding angular separation on the sky 
(a) will thus be about twice the annual parallax (since 
a = An, where A is the orbital separation expressed in 
AU). There is no indication whatsoever that the orbital 



1 The presence of Tc in HD 170970 is somewhat uncertain, 
though, since the central wavelength of the blend containing 
the A426.2 Tc line lies at the very boundary between Tc-rich 
and Tc-poor stars (Smith & Lambert 1988) 



motion of HD 35155 has been detected by HIPPARCOS. 
However, since the orbital period is of the order of the 
duration of the HIPPARCOS mission and the parallax is 
small (1.32±0.99 mas), this system represented a difficult 
challenge for the reduction consortia. The orbital motion 
has in fact been detected (flag 'O' in column 'Dupl' of Ta- 
ble |) for another short-period S star (HIC 17296 = HD 
22649; P = 596 d), but the situation is more favourable 
in this case, because of its larger parallax (tt — 6.27 mas 
compared to ir — 1.32 mas for HD 35155). 

Finally, one should be aware that the observed distri- 
bution of absolute magnitudes of a sample of stars may be 
altered by various statistical biases, depending on the se- 
lection criteria of the sample (e.g. Brown et al. 1997; Luri 
& Arenou 1997). Monte-Carlo simulations of Appendix B 
show that the observed segregation between extrinsic and 
intrinsic S stars cannot plausibly result from statistical 
biases altering the true absolute-magnitude distributions. 

4-2. Comparison with theoretical RGB and AGB evolu- 
tionary tracks 

Figure || compares the position of S stars in the HR di- 
agram with the RGB (i.e. up to the onset of core He- 
burning) for stars of different masses and of metallicity 
Y = 0.3, Z = 0.02 (Schaller et al. 1992). Figure § is the 
same as Fig. ||, but for the E-AGB, up to the first thermal 
pulse (Charbonnel et al. 1996). 

The effective temperatures of all these models have 
been converted to V — K colours using the calibration of 
Ridgway et al. (1980) for class III giants, which is strictly 
valid only in the range 2.2 < V — K < 6.8. 

The predicted location of the giant branch in the HR 
diagram is known to depend sensitively upon model pa- 
rameters like the convective mixing length or the atmo- 
spheric opacities. For the Geneva evolutionary tracks used 
here, these model parameters have been calibrated so as to 
reproduce the observed location of the red giant branches 
of more than 75 clusters (Schaller et al. 1992). The com- 
parison of these tracks with the observed location of the 
S stars in the HR diagram is thus meaningful. The main 
result of the present study is apparent on Fig. [|: the line 
marking the onset of thermal pulses matches well the limit 
between intrinsic and extrinsic S stars, so that Tc-rich in- 
trinsic S stars may be associated with thermally-pulsing 
AGB stars (the only possible exception being NQ Pup, 
but see the discussion about errors below). 

The previous result provides interesting constraints 
on the occurrence of both the s-process and the third 
dredge-up in thermally-pulsing AGB stars, by suggesting 
that those processes operate from the very first thermal 
pulses on, a conclusion already reached by several authors 
(e.g., Richer 1981; Scalo & Miller 1981; Miller & Scalo 
1982) from the luminosity distribution of carbon stars 
in th e M agellanic Clouds (see, however, the discussion of 
Sect. L3). Because there is little change in luminosity from 
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Fig. 5. The location of the RGB (up to the onset of core 
He-burning) for stars of various masses (as labelled, in Mq) 
and metallicity Y — 0.3, Z — 0.02, according to Schaller et al. 
(1992). The (V — K, spectral types) calibration is from Ridgway 
et al. (1980). Other symbols are as in Fig. ti 



one pulse to the next (see Fig. ||) , and because of the un- 
certainties affecting the location of individual S stars in 
the HR diagram, it is difficult, however, to set a limit on 
the exact number of pulses necessary to change a normal 
M giant into an (intrinsic) S star. Moreover, the statisti- 
cal biases discussed in Appendix B shift the lower bound- 
ary of the observed luminosity distribution of Tc-rich S 
stars below the true threshold. The Monte-Carlo simula- 
tions presented in Appendix B predict that the faintest 
star in the sample of 14 Tc-rich S stars displayed in Fig. [fj] 
will be observed 0.7 to 1.5 mag below the lower boundary 
of the true luminosity distribution (see Fig. Bl). This is 
well in line with the observed location of the Tc-rich S 
star NQ Pup (see Sect. Q below the TP-AGB threshold 
luminosity. 



As far as extrinsic S stars are concerned, the present 
data alone do not permit to distinguish between them pop- 
ulating the RGB or the E-AGB of low-mass stars. Nei- 
ther does the mass-transfer scenario (Sect. 1) set any con- 
straint on the current evolutionary stage of the extrinsic 
S star. However, when both the RGB and the E-AGB are 
possible, evolutionary time-scale considerations clearly fa- 
vor the RGB over the AGB. Besides, the analysis of the or- 
bital elements (Jorissen et al. 1997a) points towards them 
being low- mass stars, with an average mass of 1.6 ± 0.2 
M©. This value is in excellent agreement with their posi- 



Fig. 6. Same as Fig. but for the early AGB up to the first 
thermal pulse (Charbonnel et al. 1996). To guide the eye, a 
dashed line connects the starting point of the TP-AGB on the 
various tracks. The crosses along the right-hand axis provide 
the luminosities of the first ten thermal pulses in a 1.5 Mq, 
Z = 0.008 star (corresponding to the LMC metallicity) com- 
puted by Wagenhuber & Tuchman (1996). Other symbols are 
as in Fig. Q 



tion in the HR diagram of F ig. [s| (tw o exceptions are HD 
35155 and 57 Peg; see Sects. 4.1 and Appendix A ). 

Contrarily to what might be inferred from the smooth 
transition between intrinsic and extrinsic S stars in the HR 
diagram, the two kinds of S stars belong to distinct galac- 
tic populations, intrinsic S stars being more concentrated 
towards the galactic plane (Jorissen et al. 1993; Jorissen 
& Van Eck 1997; Van Eck & Jorissen, in preparation). Ex- 
trinsic and intrinsic S stars are thus not simply successive 
stages along the evolution of stars in the same mass range. 

Note that the lower left boundary of the region occu- 
pied by extrinsic S stars is set by the condition that T e g 
be low enough in order that ZrO bands may form. Such a 
threshold roughly corresponds to the transition between 
K and M giants (see the spectral types labelling Fig. ||), so 
that extrinsic S stars should merge into the Kill barium 
stars at lower luminosities and higher T c ft along the RGB 
(e.g., Jorissen et al. 1997a). As an example, the transition 
object HD 121447, classified either as K7IIIBa5 (Lii 1991) 
or SO (Keenan 1950), has been located in Fig. || following 
the methods presented in Sect. using the photometry 
obtained by Hakkila & McNamara (1987) and the HIP- 
PARCOS distance provided by Mennessier et al. (1997). 
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4-3. Comparison with S and carbon stars in external sys- 
tems 

Several S stars have been found in the Magellanic Clouds 
and in a few other nearby galaxies, allowing a direct com- 
parison of their luminosities with those of galactic S stars 
derived from the HIPPARCOS trigonometric parallaxes 
(Fig.0): 

• In a study of 90 long-period variables in the Magellanic 
Clouds, Wood et al. (1983) identified 14 MS stars (la- 
belled as W83 on Fig. @) in the range -7 <M hol < -5; 

• In a magnitude-limited survey of several fields in the 
outer regions of the northern LMC, Reid & Mould 
(1985) identified 10 S stars (labelled as R85 on Fig. g) 
in the range —5.0 <Mboi < —3.9; 

• In a sample extracted from the faint tail of the West- 
erlund et al. (1981) LMC survey of late-type giants, 
Lundgren (1988) finds 6 S stars (labelled as L88 on 
Fig. g) in the range -5.7 < M bo i < -4.8; 

• Bessell et al. (1983) and Lloyd Evans (1983a, 1984) 
find 16 S stars (labelled as B83, L83 and L84 on Fig. fjj 
with bolometric magnitudes provided in Tables Al and 
A2 of Westerlund et al. 1991) in LMC clusters with 
-4.8 < Mboi < -4.4 (excluding NGC 1651/3304, quite 
far from the cluster center, possibly a field star, and 
the close pair LE1+2 in NGC 1987); 

• S stars have also been found in other galaxies: for ex- 
ample, Aaronson et al. (1985) find one S star with 
M bo i = -5.15 in NGC 6822, Brewer et al. (1996) dis- 
cover one S star with Aiboi = —5.3 in M31, and Lund- 
gren (1990) finds seven S stars in the range —4.9 < 
Mbo\ < —3.1 in the Fornax dwarf elliptical galaxy (la- 
belled as L90 on Fig. g). 

The bolometric-magnitude range for all these S stars 
in external systems has been plotted in Fig. g. Their lumi- 
nosities are generally comparable to those of the galactic 
Tc-rich S stars, with the exception of the brighter W83 
S stars. These W83 S stars extend up to the theoretical 
AGB tip corresponding to the Chandrasekhar limit for the 
degenerate core. Smith et al. (1995) have shown that, in 
the Magellanic Clouds, these S stars with —7 < Mboi< —6 
are all Li- strong stars with M 4 M Q . In the solar neigh- 
bourhood, these Li-strong S stars are rare (Catchpole & 
Feast 1976; Lloyd Evans & Catchpole 1989). The only such 
star in the present HIPPARCOS sample is T Sgr, but its 
HIPPARCOS parallax is useless (Table g). 

Conversely, in all the external systems, there appears 
to be a lack of low-luminosity S stars with respect to the 
solar neighbourhood. In the Magellanic Cloud fields sur- 
veyed, this lack of low-luminosity S stars may clearly be 
attributed to the limited sensitivity of the surveys (see 
R85 and L88 in Fig. g). The situation is different in the 
Magellanic Cloud clusters, where the available surveys are 
sensitive down to Mboi = —4.0 and find many M stars but 
no S stars in the range —4.0 to —4.4. Does this mean that 
Tc-poor, extrinsic S stars are really absent from the Mag- 
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Fig. 7. Comparison of the Mboi range of S stars in the solar 
neighbourhood (this work, symbols are as in Fig. g) and in ex- 
ternal systems: 

• S stars: LMC clusters: L83 (Lloyd Evans 1983a), L84 (Lloyd 
Evans 1984), B83 (Bessell et al. 1983); Magellanic Cloud fields: 
W83 (Wood et al. 1983), R85 (Reid & Mould 1985), L88 
(Lundgren 1988); Fornax dwarf elliptical galaxy: L90 (Lund- 
gren 1990). Detection thresholds are also indicated. 

• C stars: the rightmost histogram gives the luminosity func- 
tion of the 186 LMC (thick line) and 134 SMC (thin line) 
C stars identified by Blanco et al. (1980), adopting 18.6 for 
the LMC distance modulus and 0.5 mag as the difference in 
distance moduli between the Clouds. The leftmost histogram 
provides the luminosity functions of C stars from the Wester- 
lund et al. (1991) SMC survey (hatched) as well as from the 
deeper Westerlund et al. (1995) SMC survey (open) 

ellanic Cloud globular clusters? Possibly, although extrin- 
sic S stars are mostly found at luminosities fainter than 
M bo i = —3.5 (i.e. below the RGB-tip), and were thus 
not properly surveyed in the Magellanic Cloud clusters. 
The absence of extrinsic S stars in the Magellanic Clouds 
clusters would not be surprising, though, given the sit- 
uation encountered in galactic globular clusters. Barium 
and (low-luminosity) S stars are rare in galactic globu- 
lar clusters (Vanture et al. 1994, and references therein). 
They are only present in the massive, low-concentration 
cluster cu Cen (Lloyd Evans 1983b), and their origin in 
that cluster is still unclear. Cote et al. (1997) have ar- 
gued that the binary evolution leading to the formation of 
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extrinsic heavy-element-rich stars is only possible in low- 
concentration clusters like u> Cen. In more concentrated 
clusters, hard binaries rapidly shrink to orbital separa- 
tions not large enough to accommodate an AGB star. It 
would be of interest to check whether Magellanic Cloud 
globular clusters are indeed too concentrated to allow the 
formation of extrinsic S stars. However, the recent result 
that the barium and the low-luminosity S stars in to Cen 
appear to have constant radial velocities (Mayor et al. 
1996) challenges the above picture. It may indicate that 
the barium and low-luminosity S stars in to Cen represent 
the extreme tail of the wide range in metal abundances 
observed in that particular globular cluster, as suggested 
by Lloyd Evans (1983b). Alternatively, these stars might 
have been enriched in the past by mass transfer in soft 
binary systems which were later dynamically disrupted. 

The lack of low-luminosity S stars in the Fornax dwarf 
elliptical galaxy (L90 in Fig. ^) is probably related to the 
low-metallicity of that system ([Fe/H] ~ —1.5; Lundgren 
1990). At such low metallicities, the red giant branch is 
shifted towards the blue, so that stars with luminosities 
typical of the galactic Tc-poor S stars are actually of spec- 
tral type G or K in Fornax (see Smith et al. 1996, 1997 for 
a discussion on a similar situation in the galactic halo) . No 
specific effort to find heavy-element-rich stars among the 
warm ('continuum') giants was attempted by Lundgren 
(1990). Contrarily to the situation prevailing in the Mag- 
ellanic Cloud clusters discussed above, an S star is found 
in Fornax at the luminosity threshold (M^ \ = —3.2) be- 
tween 'continuum' (G or K) giants and M giants, suggest- 
ing that the lower luminosity cutoff observed for S stars 
in Fornax is just a spectral selection effect. 

Deeper surveys are available for carbon stars, because 
of their more conspicuous spectral features. The luminos- 
ity distribution of the 134 SMC and 186 LMC carbon stars 
identified by the pioneering survey of Blanco et al. (1980) 
is shown in Fig. 0. It should be noted that the luminosity 
functions are almost identical for the LMC and SMC de- 
spite their different metallicities, so that the comparison 
with a more metal-rich galactic sample is not unreason- 
able. A GRISM survey of carbon stars in the SMC by Re- 
beirot et al. (1993; RAW) led to the discovery of 1707 such 
stars. For 100 of them, bolometric magnitudes were deter- 
mined from JHK photometry (Westerlund et al. 1991), 
with a limiting magnitude of Aiboi ~ — 3. Their distribu- 
tion is represented by the hatched histogram on Fig. fj]. 
A subsequent photometric and spectroscopic survey has 
been devoted to the ~ 5% among RAW objects having 
Mboi > -3 (Westerlund et al. 1995; Fig. g). The range of 
absolute bolometric magnitudes for carbon stars from the 
deeper SMC survey now totally covers that of galactic S 
stars, including its low-luminosity tail. 

The nature of the SMC low-luminosity carbon stars is 
still debated. They might cither be contaminating dwarf 
carbon (dC) stars from our own Galaxy, or low-luminosity 
carbon stars equivalent to the galactic R carbon stars, or 



may be extrinsic carbon stars formed by mass transfer 
across a binary system like galactic extrinsic S stars (as 
already suggested by Barnbaum & Morris 1993). Wester- 
lund et al. (1995) reject the first hypothesis, mainly be- 
cause none (except two) of their carbon stars have colours 
similar to the known galactic dC stars, and moreover, none 
exhibits a strong C2 band head at 619.1 nm, a feature ex- 
hibited by galactic dC stars. The last two hypotheses of- 
fer interesting alternatives that remain to be investigated. 
It should also be noted that the possibility that some of 
the low- luminosity SMC carbon stars be extrinsic carbon 
stars makes it difficult to use the SMC carbon-star lumi- 
nosity distribution to derive the luminosity threshold for 
the occurrence of the s-process and third dredge-up along 
the AGB (see the discussion in Sect. |4.2| ). Further observa- 
tions are clearly required in order to distinguish intrinsic 
carbon stars from possible extrinsic carbon stars. 

5. Infrared excess and position along the giant 
branches 

Infrared excesses revealing the presence of cool circumstel- 
lar material are a common feature of late- type giants (e.g., 
van der Veen & Habing 1988; Habing 1996). These ex- 
cesses are often associated with intrinsic stellar variability. 
Hacking et al. (1985) have shown that the brightest IRAS 
12/im sources outside the galactic plane are long-period 
variables (LPV). Similarly, in globular clusters, all stars 
with a 10 /im excess are LPVs (Frogel & Elias 1988), and 
moreover, they are found only above the RGB tip. These 
IR excesses are associated with a strong mass loss, which 
thus appears to be the rule among LPVs. It is expected 
to be much smaller in non-pulsating AGB and RGB stars 
(e.g., Habing 1987). Consequently, S stars with large pho- 
tometric variability and infrared excesses are expected to 
be found among the most luminous stars. Indeed, it was al- 
ready shown in Fig. || that variability in the HIPPARCOS 
Hp band tends to increase with luminosity. The present 
data offer the possibility to investigate as well the evolu- 
tion of the IR excesses probing mass loss along the giant 
branches in the HR diagram. 

The ratio R = F(12/im)/F(2.2/im), probing the pres- 
ence of cool circumstellar material emitting at 12 /im, has 
been derived from the IRAS PSC and from the K mag- 
nitude, whenever available (using the calibration of Beck- 
with et al. 1976 to convert K magnitudes into fluxes at 2.2 
fim). As already shown by Jorissen et al. (1993), there is a 
clear segregation of S stars according to their value of R. 
All Tc-deficient S stars have 0.075 < R < 0.093, consis- 
tent with photospheric blackbody colours (R — 0.073 for 
T eS = 4000 K and R = 0.093 for T cff = 3150 K, which cor- 
responds to the T c ff range spanned by extrinsic S stars). 
On the contrary, intrinsic S stars generally have R > 0.1, 
indicative of circumstellar dust. This segregation is a fur- 
ther indication of the inhomogeneous nature of the family 
of S stars. 
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Fig. 8. Infrared excess along the giant branches, as measured 
by the ratio R — F(12/j,m)/F(2.2fim). The dichotomy between 
extrinsic and intrinsic S stars is clearly apparent. Symbols are 
as in Fig. ^ 



As expected, Fig. || clearly shows that R < 0.1 indices 
are restricted to the RGB (or E-AGB), whereas S stars on 
the AGB have larger and larger R indices as they evolve 
towards cooler T e ff. 

6. Conclusions 

All stars with ZrO bands were historically assigned to a 
unique spectral class (S). However, an increasing number 
of arguments recently suggested that the S family com- 
prises stars of two different kinds: Tc-poor binary (extrin- 
sic) S stars and Tc-rich (intrinsic) S stars. The HR dia- 
gram of S stars constructed from HIPPARCOS parallaxes 
presented in this study fully confirms that dichotomy, by 
showing that: 

• extrinsic S stars are hotter and intrinsically fainter 
than intrinsic S stars; 

• extrinsic S stars are low-mass stars populating either 
the RGB or the E-AGB, and in any case are found 
below the luminosity threshold marking the onset of 
thermal pulses on the AGB. Therefore, their chemical 
peculiarities cannot originate from internal nucleosyn- 
thesis and dredge-up processes occurring along the TP- 
AGB. Their binary nature suggests instead that their 
chemical peculiarities originate from mass transfer; 

• on the contrary, Tc-rich S stars are found just above 
the TP-AGB luminosity threshold. Their location in 
the HR diagram indicates that s-process nucleosynthe- 



sis and third dredge-up must be operative quite early 
on the TP-AGB. 

The galactic latitude distributions of extrinsic and in- 
trinsic S stars are also clearly different (Jorissen et al. 
1993; Jorissen & Van Eck 1997), indicating that they are 
not simply successive stages along the evolution of stars 
in the same mass range. 
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Appendix A. 57 Peg: A very luminous extrinsic S 
star 

The star 57 Peg is peculiar in being a high-luminosity extrin- 
sic S star (M b ol= -4.3; Table §). Its outlying position can- 
not be attributed to the uncertainty on its parallax, which is 
small (ct^/tt = 0.21; Fig. §). It is a long-period binary system 
(P > 3700 d; Jorissen et al. 1997a) with a composite spectrum 
caused by a companion of spectral type approximately A3V, 
as derived by Hackos & Peery (1968) from the analysis of a 
violet optical spectrum. The A-type main sequence compan- 
ion is much better visible on the IUE spectra available in the 
archive (Peery 1986). In order to derive its exact spectral class, 
and thus its mass, its UV colours were compa red to standard 
stars from the atlas of Wu et al. (1983). Figure Al presents the 
location of 57 Peg in the ([17] - [18], [17] - [19]) colour-colour 
diagram, where [i] — [j] = — 2.5 log(F;/Fj) and F; is the average 
flux in the spectral ranges 165.5 - 175.5 nm (for i — 17), 175.5 
- 185.5 nm (i = 18) or 185.5 - 195.5 nm (i = 19). All the stars 
from the atlas of Wu et al. (1983) plotted on Fig. Al have a 



Al 



flux well above the noise level in all three bands. For 57 Peg, 
the two spectra used (SWP05384, recorded on May 28, 1979 
and SWP32554, obtained on December 18, 1987) yield iden- 
tical colours, despite the fact that they were obtained about 
8 years apart. All spectra were retrieved from the final IUE 
archive processed with the NEWSIPS software. No reddening 
correction has been applied, since it is assumed to be negligible 
given the small colour excesses listed by Wu et al. (1983) and 
by Hackos & Peery (1968). 

Main sequence stars fall along a straight line in Fig. 
with 57 Peg falling around spectral type A6. Although the 
colour-colour diagram does not probe efficiently the lumi- 
nosity class, the individual fluxes of 57 Peg are compatible 
with the companion being an A6 star on the main sequence, 
since its fluxes match those of the A6V standard HD 28527 
(SWP19459), assuming M v = +2.05 for A6V stars, and given 
the distance derived for the 57 Peg system from the HIPPAR- 
COS parallax. The corresponding mass for the companion is 
then 1.9 Mq. The S star primary must have evolved faster and 
should thus be more massive than 1.9 Mq, which is consistent 
with its position in the HR diagram along the 3 Mq E-AGB 
(Fig. ^|). That S star is thus likely to be more massive than 
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Fig. Al. The ([17] - [18], [17] - [19]) colour-colour diagram 
for 57 Peg (cross) and standard stars from the spectral atlas 
of Wu et al. (1983). Main sequence stars are represented by 
open squares with labels on the right giving the corresponding 
optical spectral type. Giants and supergiants are represented 
by filled circles, with labels on the left of the corresponding 
point 

the a vera ge extrinsic S star (for which (M) = 1.6 ± 0.2 Mq; 
Sect. |4.2| and Jorissen et al. 1997a). 

The Tc-poor nature of 57 Peg is quite puzzling, though, 
since the binary paradigm (Sect. 1) requires the companion of 
a Tc-poor S star to be a WD rather than a main sequence 
star. Possibilities to resolve this puzzle within the framework 
of the binary paradigm for extrinsic S stars include (i) 57 Peg 
is a triple system (S+A6V+WD), or (ii) the companion is an 
accreting WD mimicking a main sequence spectrum, or (iii) 57 
Peg is not an S star at all. The radial-velocity data accumulated 
so far does not allow us to test possibility (i) yet. Possibility 
(ii) is unlikely, since the SWP spectrum carries no sign of bi- 
nary interaction (like the CIV A155.0 doublet), the two spectra 
taken 8 years apart yield identical fluxes (whereas some vari- 
ability is expected in the case that the observed continuum be 
due to accretion on a WD), and finally, in the regime of rapid 
mass accretion by a WD, where a stable H-burning shell forms, 
the accreting WD mimics a supergiant rather than a main se- 
quence star (e.g., Paczyhski & Rudak 1980). Possibility (iii) - 
57 Peg is not an S star - has already been mentioned explicitly 
by Smith & Lambert (1988) with reference to an earlier paper 
(Smith & Lambert 1986), although 57 Peg is not to be found in 
that other paper. This last possibility clearly deserves a closer 
study. Otherwise, 57 Peg would add to the small set of Tc-poor 
S stars (HD 191589, HDE 332077) with a main sequence com- 
panion (Jorissen & Mayor 1992; Ake & Johnson 1992; Ake et 
al. 1994; Jorissen et al. 1997a). 



Appendix B. An evaluation of the impact of sta- 
tistical biases 

The evaluation of the statistical biases altering the true abso- 
lute magnitude distribution of any given observed sample of 
stars requires the knowledge of the selection criteria defining 
that sample. Indeed, magnitude-limited samples are mostly af- 
fected by the Malmquist bias (Malmquist 1936), whereas the 
Lutz-Kelker bias plays an important role for parallax-limited 
samples (e.g. Lutz & Kelker 1973; Lutz 1979; Hanson 1979; 
Smith 1988). For the HIPPARCOS sample of S stars consid- 
ered here, the selection criteria combine limits on the parallax 
with limits on the magnitude. Monte-Carlo simulations (see 
e.g. Pont et al. 1997) appear therefore more appropriate than 
an analytical approach to explore the biases resulting from the 
selection criteria and the observation errors when applied to a 
reasonable model of the parent population. 

The main aim of the simulations presented in this Ap- 
pendix is to demonstrate that the segregation observed in the 
HR diagram between Tc-rich and Tc-poor S stars is not an 
artefact caused by statistical biases which might in principle 
operate differentially on the two families (because e.g. of dif- 
ferent galactic distributions), but that it must result instead 
from truly different luminosity distributions. 

To this aim, a parent population of Tc-rich S stars is gen- 
erated with the following properties: 

• a large number (> 300 000) of Tc-rich S stars are created, 
with a spatial distribution characterized by an exponential 
scale height of 180 pc (as derived for the Tc-rich stars of the 
Henize sample comprising 205 S stars; Van Eck & Jorissen, 
in preparation), and a uniform projected distribution on 
the galactic plane around the Sun; 

• various bolometric-magnitude distributions are adopted 
(see below); 

• the (V — J)j colour index is derived from the bolometric 
magnitude using a fiducial AGB fitted to our sample: 

(V - 7)j = (0.5- Mbol )/l-5 + av-i gauss(0, 1) (B.l) 

where ov-i = 0.6 and gauss(0,l) is a random variable with 
a reduced normal distribution; 

• the bolometric correction is computed from the (BCi, (V — 
7)j) relation of Bessell & Wood (1984); 

• the (unreddened) Vj,o magnitude is deduced from the dis- 
tance d, Mboi, BCj and (V — J)j; 

• the Vj,o magnitude is reddened according to the cosecant 
formula (Feast et al. 1990). 

Stars are extracted from this parent population so as to 
reproduce the observed distribution of Vj magnitudes in the 
entireu HIPPARCOS sample of S stars. The "measured" par- 
allax of each extracted star is then computed from its distance 
d by 

7r mcasur od = l/d + ov(Vj) gauss(0, 1) (B.2) 

2 The V] distribution of the entire HIPPARCOS sample of 
S stars (60 stars when excluding the two close visual binaries 
and the misclassified S star; see Table |l|) has been preferred 
over the Vj distribution of the 21 S stars known to be Tc-rich, 
because of its larger statistical significance and because the 
observed Vj distributions of Tc-rich and Tc-poor S stars are 
not significantly different. 
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where the (0V,Vj) relationship is derived from a least square 
fit to the HIPPARCOS data for S stars. The simulated star 
is then retained if it satisfies the condition < o^/ti < 0.85 
imposed on the star in order to be included in the HR diagram 
of Fig. | 

In summary, the two selection criteria are: 

— Condition I: the Vj distribution of the extracted sample of 
simulated stars has to reproduce the Vj distribution of the 
HIPPARCOS sample of S stars; 

- Condition II: < g^/tt < 0.85. 

The stars extracted from the parent population and sat- 
isfying condition I (resp. I+II) define what will be called in 
the following the sample Ei (resp. Ei+u). The sample Ei is 
supposed to represent the real HIPPARCOS sample. 

Although the selection criteria that were used to include a 
given S star in the HIPPARCOS Input Catalogue are unknown 
to us, condition I ensures that they are implicitly met in our 
simulation. Besides, condition II is satisfied by 61% of the stars 
in Ei to be compared with 63% in the real HIPPARCOS sam- 
ple of S stars. This good agreement constitutes a check of the 
internal consistency of our model as a whole. 
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Fig. Bl. Bolometric magnitude distributions (upper panels) 
and their cumulative frequencies (lower panels) resulting from 
the Monte-Carlo simulations described in the text. The parent 
Mho\ distributions are rectangular and span the range indi- 
cated by the arrows in the upper panels. In the lower panel, 
their cumulative frequency is represented by the thick solid 
line. The dotted line (resp. thin solid line) refers to the dis- 
tribution of the Ei (resp. Ei+n) sample. Both histograms are 
drawn at the same scale, the dotted histogram being normal- 
ized to unity. Note that a Malmquist-type bias is responsible for 
the deviation of the dotted histogram from the original rectan- 
gular distribution, whereas a Lutz-Kelker-type bias causes the 
difference between the dotted- and the thin-solid histograms 



The comparison between the Mboi distributions of the sam- 
ple Ei+n and of the parent sample then illustrates the impact 
of the statistical biases. This is shown in Fig. Bl for rectangular 
parent distributions in the range —4 > Mboi > —6 (Fig. Bla) 
and in the range —2 > M bo i > —6 (Fig. |Bl)3). The net effect 
of these biases is to transform the original rectangular distri- 
bution into a bell-shaped distribution, with a median shifted 
by 0.4 magnitude towards fainter objets, and with asymmet- 
rical tails extending beyond the original limits. This spread is 
most pronounced on the fainter side of the distribution. The 
faintest observed S star in a sample of N objects must therefore 
be expected beyond the lower boundary of the parent distribu- 
tion (at a bolometric magnitude corresponding to the cumu- 
lative frequency 1/N). The amplitude of this offset does not 
depend on the median value of the parent distribution, but 
depends rather sensitively upon its width (compare Fig. Bl a, 



and Fig. Bl d). For reasonable values of the width of this parent 
Mboi distribution, and adopting N = 14 as for the sample of 
Tc-rich S stars plotted in the HR diagram of Fig. ^, the offset 
value is found to lie between 0.7 and 1.5 mag. According to this 
statistics, the presence of the Tc-rich S star NQ Pup at Mboi 
= —2.6 (i.e. about 1 mag below the onset of the TPAGB) is 
nevertheless fully compatible with its location on the TPAGB 
(see Fig. ||). Furthermore, Fig. Bib demonstrates ad absurdum 
that the luminosity segregation between Tc-rich and Tc-poor 
S stars is not an artefact: if Tc-rich S stars with —2 > Mboi 
> —4 typical of Tc-poor S stars were to exist, they should have 
been detected indeed. 

Two conclusions may thus be drawn from the above simu- 
lations: 

— they demonstrate that the lack of low-luminosity Tc-rich 
S stars does not result from statistical biases, since these 
biases tend to increase the number of low-luminosity stars; 

— when inferring the luminosity threshold for the appearance 
of Tc-rich S stars, it must be reminded that statistical bi- 
ases tend to make the faintest Tc-rich S stars observed by 
HIPPARCOS appear below (between 0.7 and 1.5 mag) the 
true lower-luminosity boundary of the parent population. 
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